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In a quantizing magnetic field, the chiral two-dimensional electron gas in Landau level = 0 
of bilayer graphene goes through a series of phase transitions at integer filling factors v £ [—3,3] 
when the strength of an electric field applied perpendicularly to the layers is increased. At filling 
factor zz = 3, the electron gas can described by a simple two-level system where layer and spin 
degrees of freedom are frozen. The gas then behaves as an orbital quantum Hall ferromagnet. A 
Coulomb-induced Dzyaloshinskii-Moriya term in the orbital pseudospin Hamiltonian is responsible 
for a series of transitions first to a Wigner crystal state and then to a spiral state as the electric 
field is increased. Both states have a non trivial orbital pseudospin texture. In this work, we study 
how the phase diagram at z/ = 3 is modified by an electric field applied in the plane of the layers 
and then derive several experimental signatures of the uniform and nonuniform states in the phase 
diagram. In addition to the transport gap, we study the electromagnetic absorption and the Kerr 
rotation due to the excitations of the orbital pseudospin-wave modes in the broken-symmetry states. 

PACS numbers: 73.21.-b,73.22.Gk,78.67.Wj 


I. INTRODUCTION 

In a Bernal-stacked graphene bilayer (BLG), electrons 
behave as a chiral two-dimensional Fermi gasii When 
quantized by a strong perpendicular magnetic field and 
when Coulomb interaction is considered, a rich set of 
phase transitions occurs in Landau level N = 0^^— as 
well as in higher Landau levels |7V| > 0^ The diversity 
of phases is greater in level N = 0 which has an extra 
orbital degree of freedom n = 0,1 in addition to the val¬ 
ley C = ±1 (for valleys K±) and spin cr = ±I quantum 
numbers shared by all the other Landau levels. At a fixed 
magnetic field, and in a transverse electric field E ±, a dif¬ 
ferent sequence of phase transitions is obtained at each 
filling factor iz G [—3, 3] in = 0 when a potential dif¬ 
ference Ab = E±d (or bias) between the two layers with 
separation d is increased. 

The various phases driven by the bias Ag can be de¬ 
scribed as quantum Hall ferromagnets (QHF’s). In the 
pseudospin language, the two spin, valley, and orbital in¬ 
dices are mapped into a 1/2 spin, valley pseudospin and 
orbital pseudospin respectively. In the simplest tight- 
binding model where only in-plane and inter-plane hop¬ 
pings 7o and 71 are considered (see below for a more 
precise definition of these terms) and in the absence 
of Coulomb interaction, Zeeman and bias couplings, all 
eight states in iV = 0 are degenerate. However, ad¬ 
ditional couplings such as the 74 hopping term breaks 
the orbital degeneracy while a Zeeman or a bias term 
break the spin and the valley degeneracy respectively. 
The degeneracy can also be spontaneously broken by the 
Coulomb interaction leading to different types of quan¬ 
tum Hall ferromagnetic states. At the Hartree-Fock level, 
it has been shown that Coulomb interaction completely 
lifts the degeneracy of the N = 0 octet leading to the for¬ 
mation of seven new plateaus in the Hall conductivity 
These plateaus have been detected experimentallyji 


Most of the research done so far on the QHF’s in A = 0 
has considered uniform states. But, at special filling fac¬ 
tors zz = 1,3, a sequence of phase transitions involving 
uniform and nonuniform states is also possiblei^i^ The 
nonuniform states occur in a region of bias where the 
system can be described as an orbital QHF i.e. where 
the electrons collectively condense into a linear combina¬ 
tion of the n = 0 and n = 1 orbitals. The sequence of 
transitions for zz = 3 is represented schematically in Fig. 
1 below. It involves spiral and crystal states where the 
orbital pseudospin is modulated in pseudospin space as 
well as uniform states with and without orbital coher¬ 
ence. The rotation of the pseudospin is induced by an 
effective Dzyaloshinskii-Moriya (DM) interaction due to 
the Coulomb exchange term in the Hamiltonian. There 
is no spin-orbital coupling in the 2DEC. The sequence of 
transitions in Fig. 1 is similar to that uncovered in helical 
magnets such as MnSi and Fei_a;Coa;SiiS. In these sys¬ 
tems, however, the transitions are induced by changing 
the magnetic field. 

In Refs. H and i some of us have studied several as¬ 
pects of the phase diagram in Fig. I: the order parame¬ 
ters, the density of states, the band structure, the collec¬ 
tive excitations. In the present work, we study in more 
detail possible transport and optical experimental signa¬ 
tures of the different phases in this diagram that have, 
so far, not been detected experimentally. Optical spec¬ 
troscopy has been used extensively to study graphene 
and multilayer graphene structures (see Ref. [ijJ for an 
overview of the subject). This includes absorption from 
Landau level transitions and Faraday and Kerr effects 
(optical polarization rotation of the transmitted and re¬ 
flected wave respectively). The Faraday rotation, in par¬ 
ticular, was shown to be very large in graphenei^^ Opti¬ 
cal methods can also be used to study broken-symmetry 
statesji^ Indeed, the signatures of several gapped states 
on the optical conductivity of bilayer graphene have been 
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studied in some detail in Ref. M but this did not include 
the nonuniform states that we considered in this article. 

Optical transitions between non-interacting states in 
BLG— must satisfy the selection rule \N\ —>■ |7V| ± 1 and 
are naturally classified as right and left circularly polar¬ 
ized transitions because of the change in orbital momen¬ 
tum between the two levels involved in the transition. 
The same division occurs in the simple case of = 3 in 
N = 0 where the system can effectively be mapped into a 
two-level system with orbital degrees of freedom n = 0,1 
and where spin and valley indices are polarized. (The 
phase diagram is slightly more complex for v = 1 and 
includes states with spin coherence.) In the left of the 
phase diagram in Fig. 1 (i.e. for bias < Am where 
Am defines the middle of the spiral state), electrons oc¬ 
cupy mostly the n = 0 states and optical absorption in 
all phases occurs for transition from n = 0 to n = 1. Just 
the opposite is true in the right of the phase diagram (for 
bias Ab > Am) where electrons occupy mostly the n = 1 
states. As we will show, a spontaneous uniform orbital 
QHF does not lead to absorption at finite frequency since 
its pseudospin-wave mode is gapless. The other phases, 
however, have modes that can be excited with right or 
left circularly polarized light. It should thus be possible 
to distinguish the phases by the optical absorption due 
to their collective modes. In Fig. 1, pairs of states (or 
’’conjugate states”) on each side of Am share the same 
dispersion of their collective modes as well as other phys¬ 
ical properties such as the transport gap. They absorb 
electromagnetic radiation at the same frequency although 
with different intensities and from opposite circular po¬ 
larizations of an incident electromagnetic wave. The Kerr 
rotation also show noticeably different behaviors for con¬ 
jugate biases: the polarization rotates in opposite direc¬ 
tions and the amplitude of the rotation is much larger for 
bias Ab > Am- The symmetry of the phase diagram is 
such there can be no Kerr effect at = Am- 

We extend our previous study of the phase diagraro^iSiS 
to the case where there is an electric held applied in a 
direction parallel to the layers. It was shown before that 
such a held would gap the orbital pseudospin Goldstone 
mode of the uniform phase with orbital coherencei^. We 
show in this work how the phase diagram is modihed 
by such a held, how the collective modes dispersion are 
changed in the various phases and hnally how this held 
affects the optical absorption. 

This paper is organized as follow. We describe in Sec. 
II the 2DEG at hlling factor jz = 3 as an effective two- 
level system with frozen valley and spin degrees of free¬ 
dom. We derive in Sec. Ill the phase diagram at iz = 3 
when one of the two level is hlled. In Sec. IV, we con¬ 
sider the effect of a parallel electric held on the phase 
diagram. In Secs. V and VI, we derive the optical ab¬ 
sorption and Kerr ehect due to the collective excitations. 
We conclude in Sec. VII. To avoid repetitions, we re¬ 
fer the reader to previous works in Refs. H@H, where the 
Hartree-Fock method for deriving the phase diagram and 
the calculation of the collective modes in the generalized 


random-phase approximation (GRPA) are described in 
detail. 


II. THE 2DEG AT = 3 AS A TWO-LEVEL 
SYSTEM 


The system considered in this work is a Bernal-stacked 
graphene bilayer (BLG) in a transverse magnetic B —Bz 
and electric E =Ez helds. The electric held induces a po¬ 
tential difference (or bias) Ab = Ed between the two lay¬ 
ers separated by a distance d = 3.34 A. The honeycomb 
lattice in each of these layers is described as a triangular 
Bravais lattice with a basis of two carbon atoms A„ and 
Bn, where n = 1 , 2 is the layer index and the lattice con¬ 
stant is ao = 2 . 46 A. The unit cell has four lattice sites 
denoted by {Ai, Bi, A 2 , B 2 } . The reciprocal lattice has 
an hexagonal Brillouin zone with two non-equivalent val¬ 
leys Kf = ( 1 ^) (^i’ > where ^ = ± Iji In the Bernal 

stacking arrangement, the upper A sublattice is directly 
on top of the lower B sublattice while the upper B sub¬ 
lattice is above the center of a hexagonal plaquette of the 
lower layer. 

For a neutral bilayer, the chemical potential is at the 
energy E = 0 and the low-energy excitations {E « 71 ) 
can be studied using an effective two-component model^^ 
with an Hamiltonian, in the absence of the quantizing 
magnetic field, given by 


Hln (P) = 


2m* 


Er]^P+P- 

-\aAz 


\ 

) 


( 1 ) 


This Hamiltonian is here written in the basis (A 2 , Hi) for 
valley K- and (Hi, A 2 ) for valley and p± = ± ipy. 

The parameter 


Pi 


+ 2 ^ + -), 

2 m* V 71 70 71/ 


( 2 ) 


with the effective mass m* = 2 fi^ 7 i/ 37 oaQ, where 70 = 
2.61 eV— is the in-plane nearest-neighbor hopping, 71 = 
—0.361 eV is the interlayer hopping between carbon 
atoms that are immediately above one another (i.e. 
Ai — H 2 ) and 74 = —0.138 eV is the interlayer next 
nearest-neighbor hopping term between carbons atoms 
in the same sublattice (i.e. Ai — A 2 and Hi — H 2 ). The 
parameter Jq = 0.015 eV represents the difference in the 
crystal field between sites Ai , H 2 and A 2 , Hi. We ignore 
the warping term 73, a valid approximation at the mag¬ 
netic fields considered in this articleA^d^ The Zeeman 
coupling Az — pp-bB, where 5 = 2 and <j = ±1 is the 
spin index. 

A quantizing perpendicular magnetic held is taken into 
account by making the Peierls substitution p —>■ P = 
p -f eA/c (with e > 0), where V x A =Hz. Dehning 
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the ladder operators a = {Px — iPy) £/V^ti and a'*' = 
{Px + iPy) and the magnetic length £ = yJhc/eB, 

we get 


/ 


HL = 


+ C- aa^ 
-^crAz 




C'« 

+ C+atc 


where 


c 

c± 

C' 


C ± C/3 As, 


hi 




(3) 


(4) 

(5) 

( 6 ) 


and 

hi ;* 

/3 = —= 7.24 X 10“^B[T], (7) 

71 

where the effective cyclotron frequency w* = eB/m*c 
with 


hwl = 2. 61B[T] meV. ( 8 ) 

In Eq. ([3]), the ladder operators are defined by ahn {x) = 
i\/n + l^Jn+i ix) and aipn {x) = —iy/riipn-i (x), where 
(fin (x) (with n = 0 , 1 , 2 ,...) are the eigenfunctions of the 
one-dimensional harmonic oscillator. 

The two-component model describes welli^ the 
eigenenergies and eigenstates of the Landau level N = 0 
which has eight sub-levels indexed by the quantum num¬ 
bers C and (T and an extra ’’orbital” index n = 0,1. The 
eigenenergies are 


E^,(T,n—0 — 



(9) 

- 



h— 
^ 2 

-o'^ + C + C/3Ab. 

( 10 ) 


At zero bias, the degeneracy of the octet of states in 
A = 0 is lifted by the Zeeman, the 74 and the dp terms. 
These couplings are small however since 

C = 0.39.B[T] meV, (11) 

Az = 0.12B[T] meV. (12) 

In the simplest tight-binding model where 74 = Jq = 0 
and with A^ = As = 0, the Landau le vel spectrum is 
given by En = sgn (A) y^|A| (| A| -|- l)tiuj* where A = 
0, ±1, ±2,.. .and sgn is the signum function. The gap 
between the first two Landau levels is Ei — Eq ps 37 meV. 
The two-component spinors for the levels (C, cr, n = 0,1) 
are independent of the spin index a and given, in the 
common basis {A 2 ,Bi) > by 

l!.£=+i,„,x(r) = (13) 

(14) 


with the Landau-level wave functions in the Landau 
gauge A = (0, Bx, 0) given by 

Vx (r) = (cc - A), (15) 

where X is the guiding-center index. 

Crossings between some levels of A = 0 and A = 1 
occur at bias |Ab| Ri 0.1 eV for i? = 10 T and |A_b| r; 0.2 
eV for = 30 T. Above these biases, it is in principle 
not possible to neglect Landau-level mixingj^^ 

The r = 0 K phase diagram of the chiral 2DEG in 
A = 0 has been derived in Ref. H in the two-component 
model and in the Hartree-Fock approximation for integer 
filling factors n G [—3,3] as a function of the bias Ab- 
In the present work, we study in more detail the phase 
diagram for filling factor v = S which corresponds to 
the filling of seven sub-levels in A = 0. According to 
Ref. [^, this means that all four spin up states are filled 
as well as the two states n = 0,1 with cr = —I in valley 
A_|_. Considering all these states as inert, we are left with 
a simple two-level system consisting of the two orbital 
states n = 0,1 in valley A_ with cr = — 1. We denote by 
n the filling factor of this two-level system so that n = 1 
when zc = 3. The order parameters of the different phases 
are given by 


(Pn.,n(q)) = (16) 

where is the Landau-level degeneracy and q is a two- 
dimensional vector in the plane of the 2DEG. Hereafter, 
we drop the indices ^ = —l,cr = —1 to simplify the no¬ 
tation. The diagonal elements {pn,n (q = 0)) = give 
the filling factor of each level n while the off-diagonal 
elements describe orbital coherence. 

These average electronic density n (q) is given by 


n (q) = A,^ E (-q) {Pn,m (q)) , 

n^m 

where the form factors 

3^0,0 (q) = 

A/ 1,1 (q) = e 'I 



A/ 1,0 (q) = e '’ 4 ' 
A/ 0,1 (q) = e 


({qy + iqx)i\ 

[ V2 

( Qy T ^qx) 

72 J’ 




(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


capture the orbital character of the two states. 

In a pseudospin language, the states n = 0,1 are rep¬ 
resented by the up and down pseudospin states respec¬ 
tively. The orbital pseudospin vector p (q) is related to 
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the operators Pn,m (q) by 


p(q) 

= pop (q) + Pi,i (q), 

(22) 

Px (q) 

= [po,i (q) + Pi.o (q)] /2, 

(23) 

Pv (q) 

= [po.i (q) - Pi.o (q)] /2f, 

(24) 

Pz (q) 

= [po.o (q) - pi,i (q)] /2, 

(25) 

and (p (0))| = 1/2 for z/ = 1. 



The phase diagrani at V = 1 contains both uniform 
and nonuniform states with pseudospin textures. The 
Hartree-Fock energy for a state with a pseudospin tex¬ 
ture is given by^ (apart from some unimportant constant 
terms) by 


B = 10 T and k = 5 for the host dielectric constant in 
all our numerical calculations. The different phases are 
as follows: 


0 '^O '^C A '^S* '^C* '^o* 

I Mil 



M 


Z_3 

'M 




1 

■ 


1 


L!i 


* 1* 


Ab 


FIG. 1: (Color online) Phase diagram of the 2DEG in bilayer 
graphene at filling factor P = 1. The biases Ai indicate the 
onset of each phase. (The bias axis is not to scale.) 


^ = -/3(Am-Ab)(p.(0)) (26) 

+ (P|I ■ [a (?) I + (9) A (q)] • (p|| (q)) 

^ ^ q 

+ i (^) XI (z X q) • [(p(-q)) X (p(q))] , 

X / n 


where p = {—px,Py,Pz), I is the 2x2 unit tensor and 


/ cos(2v?q) sin(2(^q) \ 

sin (2v3q) - cos {2ipq) ) ’ 


(27) 


where (/?q is the angle between the wave vector q and the 
X axis. The interactions a{q) ,b (q), c (q) and d (q) are 
defined in Ref. |^. The bias Am defines the middle of the 
spiral phase in Fig. 1. 

The interaction energy of the 2DEG with a uniform 
external electric field E|| = —VR (r) applied in the plane 
of the bilayer is given by 


(Heu) = -ey drn{r)V{r) (28) 

= - j rfrd(r) - Ell, 

with the average density n (r) defined in Eq. (II3 and 
the average total electric dipole d (q) related to the pseu¬ 
dospin operator p byi^ 

d (q) = V2e£Ny,e^^ (p|| (q)) . (29) 

To get Eq. ((28ll . we have taken into account that, for all 
phases studied in this work, the condition {p (q)) = dq^o 
is satisfied. 


III. PHASE DIAGRAM FOR v = 1 

For E|| =0, the sequence of ground states^i^i^ as the 
bias is increased is illustrated in Fig. 1. We use the values 


• Phase / for As G [0, Ao=C/i®] is a uniform 
phase where all pseudospins are oriented along z 
(i.e. i^o = 1). The Hartree-Fock gap (the en¬ 
ergy to create an electron-hole pair) is given by 
Aeh = p/2 + P (Ao — As) and is shown in Fig. 
5. Phase I has a pseudospin-wave mode with a 
dispersion loj (q) , shown in Fig. 4(a), with a gap 
oji (0) = P (Ao — Ab) /fi- This phase becomes un¬ 
stable at Ao = C/P which is the onset of the co¬ 
herent phase O. The energy 



-imeV. 

K 


(30) 


• Phase O for As G [Ao, Ao] is a uniform co¬ 
herent phase where all the pseudospins are tilted 
by an angle 6 with respect to the z axis with 
cos (9) = 1 — SP {Ab — Ao) /p- The Hartree-Fock 
gap Agfi = p/2 is constant in this phase as shown 
by the dashed line in Fig. 5. The energy of this 
phase is invariant with respect to a collective rota¬ 
tion of the pseudospins around the z axis. It thus 
has a gapless pseudospin-wave mode whose disper¬ 
sion, shown in Fig. 4(b), is highly anisotropic in 
the X — y planei^ Figure 4(b) shows the dispersion 
in the directions parallel (a;) and perpendicular {y) 
to the dipoles when they are oriented in the x direc¬ 
tion. The origin of the orbital coherence is easy to 
understand. From Eqs. (l^- (fTUl) . En=i > En=o at 
As = 0 since ^ > 0. As the bias is increased, how¬ 
ever, En=i —f En=o until En^i < En=o above a 
critical bias Ab = C/P- The electrons should then 
occupy state n = 1 instead of n = 0. But, the 
Coulomb exchange energy is smaller when the elec¬ 
trons occupy level n = 0 instead of n = 1. To re¬ 
solve this conflict, the system optimizes its energy 
by delocalizing the electrons in the two levels i.e. by 
creating a coherent state. When the small param¬ 
eters 74 , do are neglected, this coherence occurs at 
zero bias. The orbital coherence is associated with 
a finite electric polarization^^ in the x — y plane as 
shown in Eq. (Hi). Fluctuations of these dipoles 
are responsible for the electromagnetic absorption. 
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• Phase C for As G [Ac, As] is a triangular Wigner 

crystal state with one electron per site and a vortex¬ 
like texture of pseudospins around each site as il¬ 
lustrated in Fig. 2. This texture is described by the 
Fourier components (p(G)) ,where {G} is the set 
of reciprocal lattice vectors of the crystal. The crys¬ 
tal has a gapless phonon mode, which has the char¬ 
acteristic long-wavelength dispersion w ^ of 

a Wigner crystal in a magnetic field. There are also 
higher-energy gapped modes whose dispersions are 
shown in Fig. 4(c) and an electron-hole continuum 
of excitations^ The Hartree-Fock gap for this 
phase is shown in Fig. 5. 

• Phase S for Ab G [A 5 ,A 5 .] is shown in Fig. 3. 
It is a spiral phase where the pseudospins rotate in 
the plane z — n plane with n being some arbitrary 
direction in the x — y plane. The Fourier compo¬ 
nents (p (mQn)) 7 ^ 0, where m = 0, ±1, ±2,... and 
Q is the wave vector of the spiral. The energy of 
the spiral is independent of the orientation of n 
and so it has a gapless phonon mode. There are 
also higher-energy gapped mode whose dispersion 
are shown in Fig. 4(d) and an electron-hole contin¬ 
uum of excitationsi^ The dispersion of these modes 
is highly anisotropic. The gap Aet for this phase 
is shown in Fig. 5. 

The phase diagram is symmetrical with respect to the 
bias Am = C//3 + which is in the middle of the spi- 

(' 2 ') 

ral phase. We say that a phase with Ag = Am + A 

is the conjugate of that with A^^ = Am — A (with 
A > 0) in the sense that it has the same gap Aeh 
and the same spectrum of collective modes. The fill¬ 
ing factors vq > vi for A^^ and vice versa for A^^. 
We denote the conjugate phases with Ab > Am by 
In phase /*, the pseudospin are aligned 
along —z (i.e. ?'i = 1 ), the gap of the pseudospin- 

wave mode is given by w/. ( 0 ) = {Ab — Aq*) /fi and 
Aeh = b/2+/3 {Ab - Ao*) • The bias Aq. = (/P + y/Ap. 
There are no more phase transitions for Ab > Aq* in 
the two-level system. 

With the hopping parameters given in Sec. II, the 
critical biases in meV are: Aq = 53.2, Ac = 55.7, As = 
75.6, As. = 184.6, Am = 130.1, Ac- = 204.5, Aq. = 
207.0. Thus, at B = 10 T and k = 5, only the first half 
of the phase diagram falls below the limits of validity of 
the two-component model. But, the full phase diagram 
should be visible at higher magnetic field. At i? = 30 T, 
for example, Aq. < 200 meV which is within the limits 
of validity of the modekd® 


IV. EFFECT OF A PARALLEL ELECTRIC 
FIELD ON THE PHASE DIAGRAM 

Adding the dipole term of Eq. (1^ to the Hartree- 
Fock Hamiltonian and using the formalism described in 



FIG. 2: (Color online) Electronic density n (r) and dipole 
pattern (arrows) d (r) for the Wigner crystal phase at As = 
59 meV. 


n(r) 



FIG. 3: (Color online) Electronic density n(r). The arrows 
represent the dipole d (r) and pseudospin fields (pz (r)) of the 
spiral phase for As = 112 meV. 

Refs. |3@to compute the single-particle Green’s function, 
we obtain the phase diagram shown in Fig. 6. With a fi¬ 
nite parallel electric field E||, the pseudospins are pushed 
towards the x — y plane. The I and I* phases are trans¬ 
formed into the O and O* phases and orbital coherence 
is then always present. We find that the orientation of 
the wave vector Q that minimizes the Hartree-Fock en¬ 
ergy is QTE||. From Fig. 6, it can be seen that, although 
a small electric field is sufficient to suppress the crystal 
phase, a much larger field is needed to destroy the spiral 
phase. This field is 0.30 mV/nm at Ab = 0.13 eV in the 
middle of the spiral phase (not shown in the figure). The 
2DEG is described by the pseudospin energy functional 
given in Eq. (I26p where the effective pseudospin Heisen¬ 
berg exchange interaction is highly anisotropic. Figure 7 
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FIG. 4: (Color online) Dispersion relations for the (a) inco¬ 
herent phase I at As = 29 meV; (b) coherent phase O at 
As = 54 meV; (c) crystal phase G at As = 66 meV, and (d) 
spiral phase S at As = 77 meV. 

gives an idea of the strength of the parallel electric field 
which is necessary to tilt the pseudospin away from the z 
axis, against the interaction c{q) < 0 that tends to keep 
them aligned with that axis. 

When E\\ ^ 0, the Hartree-Fock electron-hole gap is 
modified in the manner shown in Fig. 5. The gap is only 
slightly increased in the incoherent and coherent phases 
and does not change noticeably in the crystal and spiral 
phases for an electric field = 0.04 mV/nm which is 
near the upper-limit of the crystal phase in Fig. 6. 

Figure 8 shows how the pseudospin-wave optical gaps 
oji (0) and wo (0) of the I and O phases change with bias 
for different values of the electric field. The electric field 
increases the optical gap in the incoherent phase I. In the 



AeCeV) 


FIG. 5: (Color online) Hartree-Fock electron-hole gap Aeh as 
a function of bias for two values of the parallel electric field 
E|| = EyX in mV/nm. The dashed lines indicate the gap in 
the uniform coherent phase O when it is not the ground state. 
This phase is replaced by the crystal and spiral phases with 
smaller gap in most of the phase diagram. 


coherent phase O, it destroys the U (1) symmetry of the 
Hamiltonian thus gapping the Goldstone modeJ^ The 
gap is practically constant in the small range of the O 
phase (the flat region of each curve in Fig. 5). Note that 
the optical gap w/ (0) —)> 0 as —)> Aq when E^^ = 0. As 
discussed below, electromagnetic absorption is expected 
at the gap frequency. 

The dispersion at small wave vector in the crystal 
phase does not change very much for an electric field 
i5|| = 0.04 mV/nm. In the spiral phase, however, the 
field can be increased to a larger value and it is possible 
to modify noticeably the dispersion. In Fig. 9, E^^ = 0 
and Aji = 0.15 mV/nm and the dispersions are shown 
in the direction of the spiral. In both the crystal and 
spiral phases, the phonon mode is not gapped by a finite 
E||. This is easily understood since the coupling to the 
external E|| involves the total electric dipole moment, a 
quantity that is not changed by a rigid translation of the 
system. 


V. ELECTROMAGNETIC ABSORPTION FROM 
COLLECTIVE MODES 

In this section, we give detailed derivation of the power 
absorbed P (w) by the collective excitations in both the 
uniform and nonuniform phases following in parts Ref. [H 
where P (w) was calculated for collective modes in quan¬ 
tum wires. The average power per unit area absorbed 
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FIG. 6: (Color online) Phase diagram of the 2DEG as a func¬ 
tion of the bias As and the electric field E||. 




FIG. 8: (Color online) Optical gaps uii (0) and ujq (0) as a 
function of bias for different values of the parallel electric field 
(in mV/nm). The flat region of each curve corresponds to the 
coherent phase O. 



FIG. 7: (Color online) In phase I, a parallel electric field Ey 
tilts the pseudospins towards the x — y plane by an angle 
0 (As, iJ||) . The field in the legend is in mV/nm. 

from an electromagnetic wave Eg (r, 2 , t) by a sample of 
size S located in the x — y plane at z = 0 is given by 

-P (w) = ^ [(E (q, z = 0, uj))* ■ j (q, w)] , (31) 

q 

where j is the surface current density in the sample and E 
is the total electric field within the sample. Since nonuni¬ 
form as well as uniform phases must be considered, it is 
necessary to use the following general relation between 
the current density and the total electric field 

j (q, w) = X ^ (q- q'. w) • E (q', z = 0 , w), ( 32 ) 

q' 


FIG. 9: (Color online) Dispersion relation of the collective 
modes in the direction of the spiral at bias Ab ~ 130 meV 
for electric field strengths = 0 and E\\ = 0.15 mV/nm. 


where a is related to the irreducible or proper part of 
the two-dimensional current response function. Equation 
(IM)l can also be written as 

j (q, w) = X ^ (q> q': w) ■ e^ (q', z = o, w), (33) 

q' 

where ^ is now related to the full (i.e. screened) 
two-dimensional current response function. In our case, 
it is the current response computed in the generalized 
random-phase approximation (GRPA). In the above re¬ 
lations, q, q^ and r are two-dimensional vectors in the 
plane of the 2DEG. The absorbed power is thus 






















where the tensor 


= (34) 

q,q' 

• ^ (q, q', w) • Ee (q', z = 0, w)]. 

From Maxwell equations, the total and external elec¬ 
tric fields E and Eg are related by 


^ (p,w) = pp-b 


^2 _ ^2p2 


(Y-pp) 


(36) 


E(q,z,a;) = Eg (q, z, w) - — y E ^ (P>‘^) (35) 

UJ Jj 

Qz 

.g.,..E'd^(q,q',a.)-Eg (q',z = 0,a;), 


with p = (q, (?z) and 1 the three-dimensional unit ten¬ 
sor. 


The total electric field is given by 


E(q,z,a;) 


Eg (q,z,w) 


47ri 1 
u) L 

47ri 1 
w L 


E 

<]z 

E 


^2 _ ^2p2 

u;2 - cV 


E^ (^’ (q'>^ = 

q' 

PP'E^ (q, q',w)-Eg (q', 
q' 


0,a;) 

z = 0, w) 


(37) 


and the absorption by 


P{^) = ^E^®[^e('4’^ = 0’‘^)’^(^’'4',w)-Eg(q',z = 0,w)] 


(38) 


q,q' 


25*2 


25*2 


E^e 


ERe 


47ri 1 




E CO 


E ^ (^’ q'> ‘^) • Eg (q', 2: = 0, w) 


q' 


47ri 1 


y_^ 


W L ^ w2 — c2p2 


E p ■ ^ (^’ (^'’ ^ 

q' 


where w is short for uj + id. 

In a crystal, translational symmetry imposes ^ (q, q',u;) —>• ^ (k-f G, k-f G', w) , where k is a vector in the 
first Brillouin zone of the reciprocal lattice and G, G' are reciprocal lattice vectors. Thus, considering a plane 
electromagnetic wave Eg (q, z = 0,a;) = SEeepSqp falling at normal incidence on the bilayer graphene system with 
polarisation vector dp, the absorption is: 


Pico) 


Ei 


Re [e; • ^ 


( 0 , 0 , ( 



47ri 1 
UJ L 

47ri 1 
UJ L 


E 

Qz 

E 


UJ^ 


— c^ql — c^G"^ 

- c^ql - c^G"^ 


\^{G,Q,uj)-ep\^ 
|g • ^ (G, 0, uj) ■ Bp 


(39) 


The lattice constant in the crystal phase is Oq = 220 A at i? = 10 T (using the relation 27rn£2 = p = 1), p follows that 
the smallest reciprocal lattice vector G « 10® m~^ and the corresponding frequency cG « 10^® rad/s. The frequency 
of the collective modes, on the other hand, is of order 10 meV i.e. « 10^® rad/s. Thus, uj << cG. Keeping terms to 
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order one in w/ cG, we get for the absorption 

~ ^R-e [e* • ^(0,0,w) • ep] 


(40) 




G^O 


27ri 

(^ V 


IcG/ 


1^ 


(G,0,a;) • Gpl-^ | V (0,0, w) • 




G 


2m 


G • ^ (G, 0, w) • Op 


The only contributions at w 7^ 0 come from the first and third terms on the right-hand side of Eq. (j40p . Defining 


^(0,0,w) = 


(41) 


we have 


Piuj^O) 


„2 77'2 „2 j ^2 

Re [d; • ^ (w) • ep] - a " 


2h 


2h 




(42) 


where a = jhc = 1/137 is the fine-structure constant and ^ (w) is now unitless. For realistic values of <5, the 
condition (uj) << 1 is satisfied and the second term in the right-hand side of Eq. can be neglected. The 

approximations that we have made to get Eq. (1^ are equivalent to neglecting retardation effects [i.e. to taking the 
limit c —^ 00 in Eq. (I39|l ]. In fact, since retardation effects are neglected in the GRPA for the conductivity, it seems 
logical to neglect them also in the calculation of the absorption. We have finally at finite frequency 


P{u;) = 


2h 


Re[e; 

r 




(w) 


(43) 


We take the incoming wave to be circularly polarized 


so that e„ 


e± 


from Eq. (H51) . 


= (x ± iy) I\p2 and the absorption is, 


P± (w) = 


2h 


Re [a± (w)] 


(44) 


where the conductivities a± are defined by 

= ^x,x P ^ ^^x,y T (^^) 

Now, the optical conductivity is related to the current 
response (t) =--E (J'J^ (r) (0)) by 


Re [cTiy (w)] = - 


Im (w)] 


OJ 


Im (w)] = 
and, for i ^ j, by 

(w) = * 


Re (w) - (0)] 


(w) 


(w -I- 

so that, for w 7^ 0 and i ^ j, we have 

R,G (^)] — : 

UJ 

Re (w)] 


Im [a,j (w)] = 


UJ 


(46) 

(47) 

(48) 

(49) 

(50) 


The total current operator is defined by J = 
—c , where is the non-interacting 

Hamiltonian of Eq. © where a Peierls substitution has 
been made to take into account an external electromag¬ 
netic field. In the two-component model, this current 
operator is related to the pseudospin operator by 


Ja ( t ) = - 


V2em^ 

ft 


(C - l3^B)Pa{T) , 


(51) 


with the convention that a = x ii a = y and vice 
versa and with pa (r) = Pa (q = 0, r). The two-particle 
current Matsubara Green’s function tensor evalu¬ 
ated at q = 0 can be related to the two-particle pseu¬ 
dospin Matsubara Green’s function tensor Ka^p (r) = 


-^{Tpa {t)pp (0)) by 


.,0{r) = -—{TJo.{T)Jp{Q)) 


(52) 






Tpa{T)pp{0) 


- ^(C-/3AB)^^,3(r). 


Finally, Ka,p ( t ) can easily be related to the two-particle 
Green’s functions defined with the operators in Eq. (HD 

i.e. 
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Xnir 


(t) = 


N 

—^ (0)) (53) 

W ni^n^) {Pn^^Ui) 


by the equations 


Kx,x = j [x++ + X+- + X-+ + X—], (54) 

Ky,y = -i [x++- X+-- X-++X—]. (55) 
Kx,y = ^ [x++ - X+- + X-+ - X—], (56) 
Ky,x = ^ [x++ + X+- - X-+ - X—], (57) 


where, to shorten the notation, we have defined the as¬ 
sociations -I- = 01 and — = 10. 

The current response is not isotropic in all phases so 
that Kx^x ^ Ky^y in general. For the conductivities at 
q= 0, 


Re [cr+ (w)] = — 
Re [a- (w)] = — 


e2 


h 


fku 


h Huj 
and for the absorption per unit area 


Im [xh— (w)] , (58) 
Im[x— h(w)], (59) 


P+(w) = - 

P- (w) = - 




2h 


fkxj 




Im [x+- (w)], (60) 
Im[x-+(w)]. (61) 


2h hw 

For linear polarizations, the absorption is given by 


For phases I and I *, the absorption is given by 


Pi') (a;) 

pP (a;) 


2h h 


P (Ao - Ab) 5 (w t w/ (0)), (65) 


e^E^.n (AB-Aof 
2h fi)^(AB-Ao.) 


i5 (w ± w/. (0))()56) 


where the frequencies w/ (0) (plotted in Fig. 8) and 
w/* (0) have been defined in Sec. III. We remark that, 
from Eq. (IMl) . it is clear that there is no contribution 
from the bubble (or polarization) diagrams in Xy± (<2-') 

so that a (0,0, w) = ^ (0,0, w) in these phases. This is 
not true in the nonuniform phases however. 

The pseudospin-wave mode in these phases is circularly 
polarized. It shows up only in P-|_ (w) for the I phase and 
only in P_ (w) for the I* phase. The response functions 

Xh _(w) and x_h (w) are equal for conjugate phases, but 

because of the prefactor in Eqs. (IBUll - dHT]) . the intensity 
of the absorption is not. The absorption P-|. (w) decreases 
as Ab —t Ao in the I phase while P_ (w) diverges (in the 
absence of disorder) at the O* —?> I* transition, decreases 
as the bias is increased and then increases linearly with 
Ab at still larger bias. 

The maximal value of the absorption frequency loj (0) 
is at Ab = 0 where hw/ (0) = 3.85 meV i.e. vi{0) = 
9.3 X 10^^ Hz in the far infrared. This frequency can be 
tuned all the way to zero by increasing Ab . It can also be 
increased by a finite Ey as shown in Fig. 8. With a finite 
E||, the absorption still appears predominantly in P_|_ (w) 
for phase I and in P_ (w) for phase /*, but the other 
circular component makes a very small contribution in 
each case. 


B. Coherent phases O and O* 


= -2h - 

X Im [x++ (w) - x+- (w) - X—t (w) + X— (w)], 


Pyi^) = - 


c^EH^-PAbY 


(63) 


2h Huj 

X Im [x-t-t (w) + X-t- (w) + X—t (w) + X— (w)], 

for an electric field polarized in the x or y direction re¬ 
spectively. 


In the coherent phases O and O *, the pseudospin mode 
is gapless and there is no absorption at finite frequency. 
As shown in Fig. 6, a finite Ey gaps that mode and makes 
it visible in absorption. The range in bias where these two 
phases are the ground state is so small that the optical 
frequency can’t change much with bias. The absorption 
is predominantly in P_|_ (w) for Ab < Am and in P_ (w) 
for Ab > Am but there is a small intensity in the other 
circular polarization which is smaller by a factor r; 10. 


A. Incoherent phases I and I* 

Analytical expressions for the absorption are possible 
in the incoherent phases when the parallel electric field 
E|| = 0. Using the results derived in Ref. H, the only 
two non-zero GRPA response functions at q = 0 in these 
phases are xh _(w) and x_h (^) with 

Xt± (^) = ^ , ■ r I T! , 1 -Y' 

Huj + i5±[C,- PAb + jWij 


C. Crystal phases C and C* 

Figure 10 shows the absorptions P± (uj) in the con¬ 
jugate crystal phases C [Fig. 10(a),(b))] and C* [Fig. 
10(c),(d)] at Ab = 66 meV and A*^ = 193.8 meV with 
A|| = 0 and A|| = 0.04 mV/nm. This figure should be 
compared with Fig. 4 (c) where the dispersion of the 
collective modes of the crystal is plotted. The gapless 
phonon mode is absent of the spectrum for w > 0. As 
in the I and I* phase, the modes active in (w) for 
Ab < Am are those active in P_ (w) for Ab > Am and 
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vice versa. The absorption peaks are more intense when 
Ab > Am however. For As < Am, the first gapped 
mode is seen only in P_|_ (w) while the fourth one is seen 
only in P_ (w) and is more than ten times smaller in in¬ 
tensity. The second and third modes are not active in 
P± (w) or Px,y (w) (not shown in the figure) for £1|| = 0. 
A finite i?|| activates the second and third modes but 
they get a very small intensity: the second mode is then 
active in P_ (w) only while the third mode, showing up 
in both polarizations, seems to be linearly polarized. The 
energy of the dominant mode is feu ~ 4 meV i.e. similar 
to the energy of the pseudospin mode at zero bias in the 
/ phase. A parallel electric field does not change notice¬ 
ably the frequency as we noted before. A linearly po¬ 
larized electromagnetic wave excites all modes of P+ (tv) 
and P_ (oj). In the crystal phases, we have checked that 
the absorption (tv) = Py (uj). 


D. Spiral phase 


Figure 1 shows the absorptions P± (tv) in the spiral 
phase at conjugate biases Ab = 77 meV [Fig. ll(a),(b)] 
and A*b = 183 meV [Fig. ll(c),(d)] with and without 
and electric field A|| = 0.04 mV/nm. This figure should 
be compared with Fig. 4 (d) where the collective modes 
of the spiral are plotted. The gapless phonon mode is 
absent of the spectrum for a; > 0. The spectrum has the 
same symmetry in P± (tv) as in the crystal phases. For 
Ab < Am, the first gapped mode is seen predominantly 
in P_|_ (w) while the almost degenerate third and fourth 
gapped mode are more active in P_ (tv). In contrast with 
the crystal phase, however, the modes of the spiral are 
not fully circularly polarized since the first (third/fourth) 
mode is only about 5 times stronger (weaker) in P+ (tv) 
than in P_ (tv). In linear polarization (not shown in Fig. 
4), both modes are stronger in Py (tv) that in P^ (y). 
Note that the spiral rotates in the z — x plane for P|| = 0 
and the absorption is maximal when the electric field of 
the electromagnetic field is perpendicular to the spiral 
i.e. in Py (tv), a fact already mentioned in Ref. The 
energy of the most intense mode is Ri 9 meV, slightly 
higher than in the crystal. At Ab = Am, the two most 
active modes are fully linearly polarized: the first in y 
and the second in x as shown in Fig. 11(e). Thus, these 
two modes become more linearly polarized as —>■ Am 

from above or from below. 

The second mode is not active in absorption when 
P|| = 0 but is activated when E\\ ^0. With a finite E\\, 
the degeneracy of the third and fourth modes is lifted 
and they become active in absorption. The second mode 
is fully linearly polarized and appears in Py (tv). 




FIG. 10: (Color online) Absorptions P± (ca) in the crystal 
phases C at As = 66 meV [(a) and (b)] and C* at the con¬ 
jugate bias As = 194 meV [(c) and (d)] with and without a 
parallel electric field = 0.04 mV/nm. 


VI. KERR ROTATION FROM COLLECTIVE 
MODES 

We adapt the calculation of Ref. [H (see Supplemen¬ 
tal material of this reference) to our specific problem to 
compute the Kerr angle for a graphene bilayer on top of 
a dielectric with a refractive index 712. We do not as¬ 
sume that (Jxx (w) = Gyy (tv) since this symmetry is not 
satisfied in all phases and do not make the simplifying 
assumptions ^axx (w) << I, (w) << 1. We take 
the incident wave E/ in the medium 1 (the vacuum) to 
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FIG. 11: (Color online) Absorptions P± (cj) in the spiral 
phases at conjugate biases Ab ~ 77 meV [(a) and (b)] and 
As = 183 meV [(c) and (d)] with and without a parallel 
electric field i?|| = 0.04 mV/nm. (e) Absorption at bias 
Am = 130 meV and for = 0. 


be linearly polarized. We write for the incident, reflected 


and transmitted (in medium 2) waves 

Ej = (67) 

Efl = + (68) 

Et = E, (t,,x + ty,y) (69) 

and for the corresponding magnetic fields 

Bj = z X E/ = (70) 


Bfl = -zxE« = E,(-r,,y + ry,x)e*(-'=i^-“*)(71) 

Bt = nsz X Et = n2Ei (4,y - ty^Sc) e*('=="-“*\(72) 

with ki = LoIc and fc2 = ^n^lc. The wave arrives at 
normal incidence on the bilayer graphene and is in part 
reflected in medium 1 and in part transmitted to medium 

2 . 

The boundary condition e| = Ef gives the equation 

E/X + El + ryjy) = E/ + tyS ), (73) 

while the boundary condition (the system is non¬ 
magnetic so that we take B = H) b| — Bf = x z 
gives 


n^Ei (txxy - tyxx) - E/y -b E/ {r^xY - Vyx^T'^) 
47r 

= —{jy^-jxY), 

where j is the induced surface current in the graphene 
bilayer. Now, according to our definition of the screened 
conductivity in Eq. dMl), we must take 

j (w) = (w) ■Ei{z = 0, w), (75) 

where ^ (w) is defined in Eq. m- We implicitly 
assumed in the above derivation that, in the crystal 
and spiral phases, we can neglect the reent components 
{j (G, w)} and keep only j (w) = j (G = 0, w) and neglect 
the diffracted components in the electric field as well. 

Solving Eqs. (TfTO) . we find for the Kerr (counter¬ 
clockwise) rotation angle from the x axis: 


tan {9k) 


Re 



2a Re 


^yx (^) 

712-1-1- ‘laOxx (w) 


(76) 


where a is the fine-structure constant. For an incident 
wave polarized along the y axis, the (clockwise) rotation 
angle is obtained by the substitution axx —f cTyy in the 
denominator of Eq. GlD. 

Figure 12 shows the Kerr angle 9k (w) for = 0 in 
the I phase and for its conjugate bias = 260 meV 
in the I* phase. We take 112 = and E|| = 0. The 
Kerr effect occurs at the frequency of the collective mode 
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which is the same for both biases but the sense of rota¬ 
tion is opposite in the two phases. The rotation in the 
I* phase takes place in a larger domain of frequencies 
(i.e. it is slower) than in the / phase and the maximum 
Kerr angle is also bigger. This is due to the prefactor 
(C — /3 Ab) in the equation for the conductivity which 
makes the conductivities much larger in phase I* and in¬ 
creases the contribution of Gxx (w) in the denominator of 
Eq. dill). 

The maximal Kerr angle is large in Fig. 12 because we 
have taken a very small value for <5 i.e. 5 = 0.014 meV. 
For 5 = 0.14 meV, the maximum angle for As = 0 is 
reduced to 4.3 degrees which is of the order of the Kerr 
angle found for other types of broken-symmetry states 
studied before.— For 5 = 1.4 meV, it is reduced to 0.7 
degrees. Since 5 approximates the effect of disorder, we 
see that the maximum Kerr angle is very sensitive to this 
parameter. 

We remark that conductivity Gxy (w) should satisfy the 
condition Re {<Jxy (w = 0)] = — 1 in the two-level system 
since v = 1. This condition is actually satisfied in phase 
I only and not in all the other phases. Indeed, Eq. (l64ll 
gives 




- 1 , 

(Ag — Aq)^ 

(As — Ao*)^ 


(77) 

(78) 


We have verified that, when levels |iV| > 0 are included in 
the calculation (in the four-component model and in the 
absence of interaction), the condition Re [axy (w = 0)] = 
—3 is then satisfied (?■ = 1 => = 3). The transitions 

from the levels N 0 occur at higher frequencies, how¬ 
ever, and they should not affect much the behavior of 
the Kerr rotation near a collective mode resonance a fi¬ 
nite frequency. 



FIG. 12: (Color online) Kerr angle for conjugate biases As = 
0 and As = 260 meV in the / and I* phases. 

Figure 13 shows the Kerr angle in the crystal and spi¬ 
ral phases for conjugate biases: (a) As = 66 meV and 
(b) Ab = 194 meV in the crystal phase and (c) A^ = 77 
meV and (d) As = 183 meV in the spiral phase. We 


take 7^2 = y/5 and = 0. As in phase I and I* discussed 
above, our calculation does not include disorder so that 
we cannot get a quantitative result for the Kerr angle. We 
use 6 ~ 0.014 meV in all curves in this figure and use the 
same step in frequency so that we can compare the rela¬ 
tive size of the rotation angle. As in Fig. 12, we expect 
that if disorder is included, it will reduce the Kerr angle. 
The modes visible in the Kerr rotation are also those ac¬ 
tive in absorption as can be seen by comparing Fig. 13 
with Figs. 10-11. In the crystal and spiral phases where 
two modes are active in the range of frequency shown in 
Fig. 12, the Kerr angle for Ab < Am increases in going 
towards the resonance in the lower-energy mode while it 
decreases when approaching the fourth mode (the second 
and third mode do not appear in the Kerr rotation). The 
rotations are in the opposite directions for the conjugate 
biases. The Kerr rotation is bigger by a factor Ri 10 in 
the conjugate phases with Ab > Am- Since the same 
mode leads to opposite rotation of the polarization in for 
conjugate biases, the rotation must cease at bias Am . We 
have checked numerically that this is the case. 

In all phases but the spiral phase, we get numerically 
that the Kerr angle is the same for an incident wave po¬ 
larized along the x or y axis. This is also true at small 
bias {Ab < Am) in the spiral phase: a y polarization 
gives the same result as in Fig. 13 (c). For Ab > Am, 
however, the rotation angle is bigger when the polariza¬ 
tion is along the x axis (i.e. for a spiral rotating in the 
z — x plane) as can be seen by comparing Fig. 13 (d) and 
Fig. 13 (e). 


VII. CONCLUSION 

We have continued in this work an analysis of a 
sequence of phase transitions involving uniform and 
nonuniform states with orbital coherence initiated in Ref. 
1^. This sequence, represented in Fig. 1, occurs at filling 
factor V = Z oi bilayer graphene, in Landau level A = 0, 
when an electric bias between the two levels is increased. 
At filling factor zz = 1, the same states are found in the 
phase diagram but others involving spin coherence are 
present as well. In this work, we choose to concentrate 
on the simpler case of iz = 3 in order to simplify the 
discussion but the same results would apply to the cor¬ 
responding states at v = 1. 

Because we were only interested in studying the sig¬ 
nature of the collective excitations due to transitions be¬ 
tween the two levels n = 0,1 in Landau level A = 0, we 
could use the two-component modeli^ to simplify the cal¬ 
culations instead of working with the full four band model 
of bilayer graphene. All calculations were done at T = 0 
K using the Hartree-Fock approximation to derive the 
phase diagram of the interacting chiral electron gas and 
the generalized random-phase approximation (GRPA) to 
obtain the dispersion relations of the collective modes. 

Our calculations show that there are qualitative differ¬ 
ences in the transport gap and in the optical properties 







14 


(b) 


10.0 - 
5.0 


(C) 1.0 

1 »- 0-5 

a> 

o 

O) 

“ 0.0 
C 

- 0.5 

- 1 . 0 , 


(d) 


20 


05-10 



- 

\ Ag=194meV 






V 

■ \ 



\, 




6 8 
(0 (in meV) 


10 


J 

Ab= 77 meV 




10 12 
10 (in meV) 


- 

Ag=183 meV 

X polarization 

- 



10 12 
(0 (in meV) 


16 



(electromagnetic absorption and Kerr rotation) in the 
different phases represented in Fig. 1. The absorption 
frequency can be tuned all the way to zero in the inco¬ 
herent phase while it is almost constant in the crystal and 
spiral phase and zero in the uniform state with orbital co¬ 
herence. The absorption frequency can be modified by 
applying an electric field in the plane of the layers, but 
in a limited way since such a field suppresses the nonuni¬ 
form states. Another clear signature occurs in what we 
call the conjugate states where the same optically active 
modes at two conjugate biases are active in opposite cir¬ 
cular polarizations. Moreover, the optically active modes 
in the incoherent and crystal phases are circularly polar¬ 
ized in contrast to the active modes in the spiral phase 
which are neither completely linearly polarized nor com¬ 
pletely circularly polarized. 

Another qualitative difference that we find is in the 
behavior of the Kerr rotation. There is a Kerr effect 
near the frequency of each mode active in the absorp¬ 
tion. Conjugate states show a Kerr effect at the same 
frequency but with opposite sign for the rotation of the 
polarization. The maximal polarization angle (near the 
resonant frequency) is bigger for bias A b above the mid¬ 
dle of the spiral phase. The observability of the Kerr 
effect would depend very much on the strength of the 
disorder in the sample. A proper treatment of disorder is 
also necessary to derive the actual absorption line shape 
due to the collective excitations. 

The sequence of phase transitions that we studied in 
this paper has not been observed so far. We think that 
the work presented here may help in distinguishing the 
different phases. Of course, several improvements would 
be necessary to get numerically accurate predictions. 
Landau level mixing, for example, is known to modify the 
dispersion and hence the absorption frequency of the col¬ 
lective modes or magnetoexcitons^S. Also, an important 
information as regards the experimental observability of 
the phases would be to know how fragile they are with 
respect to thermal fluctuations. We leave these problems 
for further work. 
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